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Abstract. We consider multiple orthogonal polynomials with respect to two 
modified Jacobi weights on touching intervals [a, 0] and [0,1], with a < 0, 
and study a transition that occurs at a — — 1. The transition is studied in 
a double scaling limit, where we let the degree n of the polynomial tend to 
infinity while the parameter a tends to — 1 at a rate of 0(n _1//2 ). We obtain 
a Mehler-Heine type asymptotic formula for the polynomials in this regime. 
The method used to analyze the problem is the steepest descent technique for 
Riemann-Hilbert problems. A key point in the analysis is the construction of 
a new local parametrix. 

1. Introduction and statement of results 
1.1. Introduction 

Multiple orthogonal polynomials are a generalization of orthogonal polynomials 
that originated in works on Hermite-Pade rational approximation problems, but 
recently found other applications in random matrix theory and related probabilistic 
models. 

In the approximation theory literature two main classes of multiple orthogo- 
nal polynomials were identified for which detailed asymptotic results are available. 
These are the Angelesco systems and the Nikishin systems. In an Angelesco system 
[I] the multiple orthogonality is defined on disjoint intervals, while in a Nikishin 
system [3 2) the orthogonality is on the same interval with orthogonality measures 
that are related to each other via an intricate hierarchical structure. 

A main stimulus for the asymptotic analysis of orthogonal polynomials was 
given by the formulation of a 2 x 2 matrix valued Riemann-Hilbert problem for 
orthogonal polynomials by Fokas, Its and Kitaev [21] and the subsequent applica- 
tion of the powerful Deift-Zhou steepest descent technique to this Riemann-Hilbert 
problem in [15l [16] and many later papers. 

A Riemann-Hilbert problem for multiple orthogonal polynomials was formu- 
lated by Van Assche, Geronimo and Kuijlaars 41 . The Riemann-Hilbert problem 
is of size (r + 1) x (r + 1), where r is the number of orthogonality weights for the 
multiple orthogonal polynomials. The Riemann-Hilbert formulation was already 
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used in several papers, see e.g. [U EJ El EH HU H2 E3 H3 EI] for the asymp- 
totic analysis of multiple orthogonal polynomials and their associated multiple 
orthogonal polynomial ensembles [SSI ES] ■ 

In this paper we consider Angelesco systems on two touching intervals [a, 0] 
and [0,1] with a < —1. Our interest is in the special behavior at that takes 
place near a critical value of a. A prime example for this situation is given by the 
Jacobi- Angelesco weights 

w 1 {x) = \x-a\ a \x\ f3 \x-l\\ x£(a,0), 

w 2 (x) = \x-a\ a \xf\x- z 6(0,1), lJl! 

with a, /3, 7 > —1, which were hrst studied by Kaliaguine [331 HI]- The associ- 
ated multiple orthogonal polynomials are among the classical multiple orthogonal 
polynomials [5] and as such have a number of very special properties. There is 
e.g. a raising operator which gives rise to a Rodrigues-type formula and a third 
order linear differential equation as well as an explicit four term recurrence rela- 
tion for the diagonal case Jacobi-Angelesco multiple orthogonal polynomials, see 

[21 12H EH EH E3 EH SO] - 

1.2. Modified Jacobi-Angelesco weights 

We generalize the system (ll.lj) by considering more general modified Jacobi weights 
on the two intervals (a, 0) and (0, 1). We will use the following weights w\ and w 2 
throughout this paper. 

Definition 1.1. Let a < 0, a, (3, 7 > — 1 and define 

A 1 = [a,0], A 2 = [0,1]. 

For j — 1,2, let hj be strictly positive on Aj with an analytic continuation to a 
neighborhood of Aj in the complex plane. Then we define 

wi(x) = (x - a^lxfhxix), xGAi, 

w 2 (x) =x?{l-xyh 2 (x) 7 xeA 2 . (1 ' 2) 

When appropriate we set Wj (x) = for x G E \ Aj . 

The definition of the multiple orthogonal polynomial (of type II) with respect 
to the weights (|1.2j) is as follows. 

Definition 1.2. Given a multi-index (ni,n 2 ) £ N 2 the multiple orthogonal poly- 
nomial is defined as the unique monic polynomial P ni .n 2 °f degree n% + n 2 such 
that 

/ ^ni,n 2 {x)x Wj(x) dx = 0, for k = 0, . . . rij — 1, (1-3) 

JA-j 

for j = 1,2. 

Since we are dealing with an Angelesco system of weights [T| the polynomial 
P ni ,n 2 indeed exists and is uniquely characterized by (|1.3|) . It is also known that 
all the zeros of P ni ,n 2 are rea l an d simple with ri\ zeros in (a,0) and n 2 zeros in 
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(0,1), see e.g. [39]. For the definition of the multiple orthogonal polynomials of 
type I we also refer to [39] . 

1.3. The phase transition 

We consider in this paper the diagonal case 

m = n 2 = n. 

It is known that the zeros of the multiple orthogonal polynomial P n , n have a weak 
limit as n — » oo, which only depends on the parameter a < 0. The limiting zero 
distribution can be characterized as the solution to a vector equilibrium problem 
for two measures [22J, [32] ■ 

Define the logarithmic energy I(u) of a measure v as 

I{v) := [f log ^— dv(x)dv(y), (1.4) 
J J \x-y\ 

and the mutual logarithmic energy I(u, fi) of two measures v and [i as 

I{u,fi) := log dv(x)dn{y). (1.5) 

Then the vector equilibrium problem is defined as follows. 

Definition 1.3. The vector equilibrium problem asks to minimize the energy func- 
tional 

E(v u v 2 ) := J{vi) + I{vi,u 2 ) + I(y 2 ) (1.6) 
among positive measures v\ and v 2 with supp(^i) C [a,0], supp(^2) C [0,1] and 
fdut = i, $dv 2 = \ 

One may interpret this energy functional as the energy resulting from two 
conductors [a, 0] and [0, 1] with each an equal amount of charged particles. Parti- 
cles on the same conductor repel each other, such that the resulting electrostatic 
force is proportional to the inverse of the distance between the two particles, which 
accounts for the terms I(vi) and I(v 2 ). Additionally, particles on different conduc- 
tors also repel each other, but with only half the strength. This leads to the term 

v 2 ) in (|1.6p . This kind of interaction is known as Angelesco-type interaction, 
see [3]. 

The minimizers v\ and v 2 for the Angelesco equilibrium problem are called the 
equilibrium measures. They exist, are unique and are absolutely continuous with 
respect to the Lebesgue measure. It is due to Kaliaguine [23 that the endpoints 
of the supports of the equilibrium measures are given by a, 0, 1 and a fourth point 
b: 

(a + lf 



9(a 2 -a + 1) 

such that 

supp^i = [a,b] C [o, 0], suppz^ = [0,1], if a < —1, 

supp^i = [a, 0], supp^2 = [6,1] C [0,1], if a > — 1. 



(1.7) 
(1.8) 
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Figure 1 . A sketch of the equilibrium densities in the three cases 
a < — 1, a = — 1, and —1 < a < 0. For a clearer picture the size 
of the gap between and b has been exaggerated. In a true plot 
b would be much closer to as the gap is less than one ninth of 
the length of the larger interval. 

We see here the pushing effect: the charge on the smaller interval pushes away 
the charge on the larger interval, thereby creating a gap in the support. The gap 
disappears in the symmetric case a = — 1 where we have 6 = 0. 

The density of the equilibrium measures blows up as an inverse square root 
at the endpoints a, 0, and 1 of its supports. These are the so-called hard edges. 
For o -1 there is a soft edge at b, where the equilibrium density vanishes like 
a square root. In the symmetric case a = — 1, where both intervals have equal 
size, both measures have full supports and the densities behave like inverse cube 
roots at 0. We may call in this case a Kaliaguine point after [23]. A sketch of the 
densities in the three cases is given in Figure [1] 

It is this interior soft-to-hard edge transition as a varies around — 1 that will 
give rise to a new critical behavior of the multiple orthogonal polynomials around 
that we wish to describe in this paper. 

1.4. Main result 

The main result of this paper is a Mehler-Heine type asymptotic formula for the 
multiple orthogonal polynomial P n ^ n (z) near z = 0, with the parameter a near — 1. 
We use P n ^ n (z; a) to denote the dependence on a. 

Theorem 1.4. For a < close enough to —I, let P n ^ n (z;a) be the multiple orthog- 
onal polynomial with respect to the weights (|1.2|) and the multi-index (n,n). Let 
t € M and _ 

a n = -l + ^-. (1.9) 

Then we have for every z£C, 

where C n is a positive constant and 

Q(z;T)=i^t-P- x exp(-^-^-+?J dt, zeC, (1.11) 



Double scaling limit for modified Jacobi-Angelesco polynomials 



5 




Figure 2. The contour Tq appearing in the Mehler- Heine for- 
mula (|1.10[) . The dashed line denotes the cut of t^P^ 1 . 



where the contour Tq is shown in of Figure® The convergence in (|1.10|) is uniform 
for z in compact subsets of C. 

The function Q in is an entire solution of the third order differential 

equation 

z 2 Q"'(z) + 2(/3 + l)zQ"(z) + ((3 2 + /3 - tz)Q'(z) + (z - rp)Q(z) = (1.12) 

The differential equation has a regular singular point at z = 0, with associated 
Frobenius indices equal to 0, — (3 and —,0+1. For (3 > — 1, there is a one-dimensional 
space of entire solutions to (|1.12l) . unless (3 = in which case this space is two- 
dimensional. In case (3 ^ we may characterize Q as the unique entire solution of 
(|1.12j) satisfying 

0(0) = j^. (1-13) 
The constant C n in (|1.10[) is given by 

c "=wf5^- tV+v *" 4 (^)" (li4) 

where c\ and C2 are positive constants defined in (|4.44|) below. They are determined 
by the analytic factors h\ and hi in the weights (|1.2|) . and are independent of n. 
For simple analytic factors hi and hi one can evaluate c\ and C2 explicitly. For 
example, if hk is a constant function, then Ck = 0. 

Remark 1.5. In the case r = the function Q(z;t) from (jl.lip can be written as 
a generalized hypergeometric function 

n( n\ 2tt / (3 + 1 (3 + 2 z 2 

Q( z i 0) = 7T71TTTT 0^2 



r(/? + i) 

This function was already found by Sorokin [34] in a Mehler-Heine formula for 
certain multiple orthogonal polynomials of Laguerre-type. More recently, it was 
obtained for Jacobi-Angelesco multiple orthogonal polynomials by Tulyakov [37] 
and Takata [3B], who both prove Theorem II .41 for the case r = and weights (|1.2[) 
with hi = 1, /12 = !■ 
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Remark 1.6. In [40] an explicit formula for the Jacobi-Angelesco polynomial P^rJ 3 ' 1 ^ 
is given, that is, the multiple orthogonal polynomial with weights (|1.1[) . namely 

( 3n+ °„ + " +7 )p&"-W) 

Applying Stirling's approximation formula to the binomial coefficients, one can 
then derive that with a n given by (|1.9[) . 



P„,»(0;a n ) = (-l) r 



2tt 1 /2^ Q+ ^+ 7 



r(l + )8) V&F V3 



-♦(£)"(. + o (.-»)) ,00, 



which is consistent with (|1.14j) . since in this case one can evaluate ci and c-i to be 
Cl = 7 (log2-±log3), c 2 =a(log2-±lo g 3). (1.17) 
1.5. Overview of the rest of the paper 

We use two main tools to prove the Theorem ll.4| namely Riemann-Hilbert (RH) 
problems and modified equilibrium problems. These will be discussed in the next 
two sections. 

The RH problems are of size 3x3. We first discuss the RH problem for 
the multiple orthogonal polynomial P ni ,n 2 with the modified Jacobi weights. The 
steepest descent analysis of the paper will lead to a local parametrix that is built 
out of a local model RH problem, that is discussed in detail in Section |2~2"1 This 
model RH problem is new, although it is related to another model RH problem 
studied recently in a different connection 28 . 

The modified equilibrium problem is related to a Riemann surface in Section 
[3] The same Riemann surface will also play a role in the construction of the outer 
parametrix in the steepest descent analysis. 

Section [4] is the bulk of the paper. It contains the steepest descent analysis 
of the RH problem for multiple orthogonal polynomials. It follows the usual steps 
in such an analysis as e.g. done in O [X5l [29]. In a first transformation we use 
the g-functions coming from the modified equilibrium problem to normalize the 
RH problem at infinity. The next transformation is the opening of lenses. Then we 
construct outer and local parametrices that are used in the next transformation. 
After this transformation one typically arrives at a RH problem that is normalized 
at infinity, and for which the jump matrices all tend to the identity matrix as 
n — > oo. It is a curious fact that this does not happen in the present paper. The 
jump matrix on a circle around will take the form 

I + Z n (z) + 0(n-V 6 ) 
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where Z n (z) is bounded on the circle, but it does not tend to as n — > oo. We can 
resolve this problem by making another transformation, where we use the special 
structure of the matrices Z n (z). This extra step in the steepest descent analysis is 
also used in the recent papers []j|J [2E] , which makes it reasonable to suspect that 
the need for such an extra step is a more common phenomenon in the steepest 
descent analysis of larger size RH problems in a critical situation. 

The proof of Theorem 11.41 is given in the final section |6l Here we unravel all 
the previous transformations, and we pay special attention to the behavior around 
0. 

In a forthcoming paper we plan to analyze the determinantal point process 
that is associated with the modified Jacobi-Angelesco weights. This is an example 
of a multiple orthogonal polynomial ensemble |25] where half of the particles are 
on [a, 0] and the other half are on [0, 1]. There is again a critical behavior at as 
a varies around — 1, and we will find a new family of limiting correlation kernels 
in this setting that are also related to the solution of the local model RH problem. 



2. First tool: RH problems 
2.1. The Riemann-Hilbert problem 

The multiple orthogonal polynomial (|1.3|) are characterized in terms of a 3 x 3 
matrix valued Riemann-Hilbert problem (RH problem) due to [41]. We use the 
RH problem for the asymptotic analysis to derive our results. 

We work with the modified Jacobi weights w\ and u>2 (| 1 . 2|) and we take 
a general multi-index (711,712)- The RH problem then asks for a function Y : 
C\ [a,l] -> C 3x3 such that 

• Y is analytic on C \ [a, 1], 

• Y has continuous boundary values Y± on (a,0) and (0, 1) satisfying a jump 
relation Y + = Y^ Jy with jump matrix 

wi(x) w 2 (x)^ 

J Y (x) = I 1 ), i€ (a,0)U (0, 1), (2.1) 

1 

where it is understood that wi(x) = on (0, 1) and W2(x) e on (—a, 0), 
Y has the asymptotic behavior 



(2.2) 




as z 00, and 
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• Y has the following behavior at the endpoints of the intervals 
Y(z) = C ' | <(:) | as ; 




where e(z) 



< 



\z-a\ a ifa<0, 
bg|z-a| ifa = 0, (2.3) 
1 if a > 



Y(z) 




where e(z) 




7<0, 
7 = 0, 
7 > 



(2.4) 



e(z) e(*)> 
Y(z) = 0\ l e(z) e(z) | , as z -> 0, 
e(z) e(z). 



where e(z) 




if 

if = 0, (2.5) 
if p > 0, 



where the O is taken entry-wise. 

As in [TH |41] one can show that there is a unique solution of the RH problem, 
see also [29] for the role of the endpoint conditions (|2.3[) . (|2.4[) and f|2 . 5[) . The first 
column contains the multiple orthogonal polynomials of type II with respect to 
multi-indices (m, 112), (n\ — 1, TI2) and (ni, ri2 — 1) and the other columns contain 
Cauchy transforms of the polynomials times the weights. Indeed, the solution is 
equal to 



P 



ill ,n 2 



1 f0 ^i.^W^W ^ 
27TZ J a x — z 

Y(z) = I d lPni ^ n2 (z) £-J° a ^-W*W*) 



K <kP ni ,n 2 -l(z) 2^l! a 



P„ 



x — z 



dx 



,-i(x)wi(x) 



1 r J 

27TI JO 

H JO 



1 P n 



(x)w 2 (x) 



1 -P, 



277; 



^ 2m /o 



1 -P, 



x — z 

■■1-1, n 2 Wg2W 

X — Z 
■ 1 ,n 2 -l(x)w 2 (x) 



da- 



da: 
dx j 
(2.6) 

for certain non-zero constants di and c?2- The inverse matrix Y 1 contains multiple 
orthogonal polynomials of type I. 

We apply the Dcift-Zhou steepest descent analysis to the RH problem for Y 



with n\ — ri2 = n in the limit where n — ¥ 00 and a 
a number of transformations 



-1 



v/2t 



-1. Via 



Y h^T^ R 
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Figure 3. The contour E$ and the jump matrices of ty. 

we arrive at a matrix valued function R that tends to the identity matrix as 
n — > oo. 

Particularly relevant references on the steepest descent method for this paper 
are [TU1 HS1 [55] , see also [TT]. 

2.2. The local model RH problem 

At a crucial step in the steepest descent analysis we need to do a local analysis at 
the point 0. We have to construct there a local parametrix that will be built out 
of certain special functions. In non-critical situations this can be done with Bessel 
functions of order /3, but in the critical regime that we are interested in we need 
functions that satisfy a third order linear differential equation. They are combined 
in a 3 x 3 matrix valued RH problem that we call the local model RH problem 
and that we describe next. 

In the local model RH problem we are looking for a 3 x 3 matrix valued 
function that depends on two parameters /3 > — 1 and r G R. Since (3 is con- 
sidered fixed, we do not emphasize the dependence on (3. We may write ^(z; r) to 
emphasize the dependence on r. Then should satisfy the following. 

• VP is defined and analytic in C \ where is the contour consisting of the 
six oriented rays through the origin as shown in Figure [3l 

• \& has continuous boundary values on that satisfy the jump condition 

where the jump matrices J* are also given in Figure [3J 
Thus the parameter j3 appears in the jump condition. 

The dependence on r is in the asymptotic condition as z — ¥ oo. We write 
throughout this paper 
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r,3 o 



ri 




Figure 4. The contours Ti, T 2 and T 3 in the t-plane. The dashed 
line denotes the cut of t~@~ 3 . 
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• As z — > oo with ± Imz > 0, we have 

= sj^e^z^ 1 ) n± (/ + O B±e e ^, (2.7) 

\ z~ s/ 

where f2±, £?± and 6(2; t) are defined by 

/-w 2 1 lj\ 
0+ := I 1 — 1 —II, B + := 

\ n n p-^f 1 / 

(2- 

/e~T- \ 

<).:= I - I -I —J I . B- : = 

and 

9( Z . T ) •= fdiag (^1(2; r), e 3 (z; r), 6*2(2; r)) for Imz > 0, 
\diag (6*2(2;; t^^-t^i^t)) forlmz<0, 

and the 6^ are defined by 

6> fe (z;r) := ~uj k z% - ruj 2k z^ forfc = l,2,3. (2.10) 

The expansion (|2.7j) for ^(z) as 2 — > 00 is valid uniformly for r in a bounded 
set. 

We construct ^(z; r) out of solutions of the third order linear differential 
equation 

z q i"(z) - pq"{z) - rq'(z) + q(z) = 0. (2.11) 

Note that this is not the same differential equation as (|1.12[) . However, the two 
are related, since if q satisfies (12.111) then 

Q(z) = z-f j q"(z) 

satisfies (1TTT2")) . 
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The differential equation ()2.11[) has solutions in the form of contour integrals 



:--V+zi 



(2.12) 



where T is an appropriate contour so that the integrand vanishes at the endpoints 
of the contour T. Define three contours ri,r 2 and T 3 as in Figure HI and define 
for z with Re z > 



qj(z) = / t 



f-/3-3. 



r+zt 



dt, 



3 = 1,2,3, 



(2.13) 



where we choose the branch of t " 3 with a cut on the positive real axis, i.e., 



-,3-3 _ \ t t-P-3 e (-P-3)iaxgt^ 



< arg t < 2tt. 



The integrals (j2.13[) only converge for z with Re z > 0, but the functions qj can 
be continued analytically using contour deformations. Branch points for the qj- 
functions are and oo and we take the analytic continuation to C \ (-co, 0], thus 
with a branch cut on the negative real axis. 



Definition 2.1. Define W in the upper half plane by 





(e^qt 


e^ l q 3 qi\ 






e 2 ^q' x 


q> 2 






\e 2 ' 3m q'{ 


e^q's' <&) 






f e ^i qi - 


f 92 e^q 3 


92 \ 




e 2 ^q'i - 


f q' 2 e^q' 3 


92 , 




\#P*i<£ - 




9 2 7 






f q 2 - e 2 ^ l q z 


e?™q 3 




e 2 P™q'i ' 


f q' 2 - e 2 ^ 


e^ l q' z 




\e^ m q'{ - 


f q" _ e ^ q l 


e 0ni q n 



< arg z < 



4 - 



92 
92 



§ < argz < ^ 



2f < argz < 7T, 



(2.14) 



and in the lower half plane by 





( 92 


e /37r 


93 














92 


e /37T 


9.3 


-e 2 ^ 








-f < argz < 0, 








93' 


-e^q'l) 












^92 - 






e^q 3 


-e 2/37r 










9 2 - 




rig/ 




_ e 2/37r 






< argz < -f 




U- 




9i 


e /3ri g // 


_ e 20m 


9iV 










"91 - 


f 92 


+ 93 e^q 3 - 


e 2/3 


™9i^ 






e 2/3 


™9i - 


f 92 


+ 9 3 e^Y 3 - 


e 2/3 


ri 9 / 


, -7r<argz<-^ 






f 92 


+ q 3 ' 


93 


e 2 ^ 







(2.15) 



It is then an easy exercise to check that ^> indeed satisfies the required jumps 
= "J- J* on the rays argz = 0, ±j,±^-. For the jump on the negative real 
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axis however, we have to take into consideration the behavior of the functions qj(z) 
as z circles around 0. Using contour deformations one can show that for z < 0: 

qi , + (z)\ /l + e 2 ^ 1 0\ 

&,+ (*) = -e 2 ^ q 2 ,-(z) . (2.16) 

<&,+(*)/ V <? Pni i V W(*V 

The jump of $ on the negative real axis follows from this in a straightforward way. 
As for the asymptotic behavior, we have 

Proposition 2.2. The function \& defined in (|2.14l) and (|2.15[) satisfies the asymp- 
totic condition (|2.7p . 

Proof. This follows from a classical steepest descent analysis applied to the contour 
integral representations for the qj (|2 . 13[) . Define the phase function z, r) by 

0(t;z, T ) : =I--L + zt (2.17) 

The main contributions in the integrals occurs around the saddles t = t k = 
tfc(z; r) of which are the solutions to 9'{t) = 0: 

t k (z;r) = -Lu k z-? - ^u 2k z-i +0 (z - ') , fc = l,2,3. (2.18) 

The O-term here is uniform for r in compacta. The critical values are given by 

9{t k ) = 6(t k ;z, r) = ~J : z§ - ™ 2 M + L. + O . (2.19) 

We also need the second derivative of 6> in the saddle points: 

9"(t k ) = -3w 2 M + O(z). (2.20) 

Through each saddle point t k there is a steepest descent path Tfe. This is a 
path such that lm9(t) — Im9(t k ) for all t G Y^. Let o^, |aj.| = 1 be the tangent 
direction of T k in t k . The steepest descent method then yields 

f t-"-VW d£ = a J _„J^ .. 2 t k (z)-^e s ^ (l + ©(z- 1 ^)) . 

(2.21) 

The fact that is a steepest descent path guarantees that 9"(t k (z))al is 
negative. Substituting ([2~T5]) . ([2"TH?J) and (j2~2TJ)) we find 

= ±Me*uPz-H-utz-i)-'-*e-i" h *''-™ ak > 1 '' (l + O (*"*)) . 

The final step in the proof is the identification of the steepest descent paths, 
and the deformation of the T k into steepest descent paths. This gives us expressions 
for the q k in each sector. For the second and third row of '3/ we remark that by 
()2.13p differentiation of the q k is equivalent to increasing /3 by 1. 
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The final expansion for ^ then turns out to be exactly as in ()2.7[) . □ 

For the further analysis we also need to know the next order term in the 
expansion (|2.7[) . 




Lemma 2.3. We have as z — > oo with ± Imz > 
tf(z) 

V ° \0 

x (l+ (*i)±*-* + O (*"*)) B ± e e(z \ (2.22) 
where the constant matrices (^i) ± are given by 



(2.23) 




(2.24) 



Proof. The first row of \&i can be found in a straightforward way by expressing 
that the entries in the first row of &(z) solve the differential equation (|2.11l) . For 
the second and third rows we increase /3 in (12. lip by 1 and 2, respectively. □ 



3. Second tool: modified equilibrium problem and Riemann surface 
3.1. Modified equilibrium problem 

One of the transformations in the Deift/Zhou steepest descent analysis of the RH 
problem for (multiple) orthogonal polynomials is typically based on the limiting 
zero distribution of the associated polynomials. In our situation these are given 
by the vector equilibrium problem from Definition II. 31 As explained in Subsection 
11.31 the endpoint b = b a of one of the supports is varying with a and tends to 
as a — > — 1. Working with measures with varying supports around would 
cause major technical problems. Therefore, following |1Q[ 112 ) . we use a modified 
equilibrium problem, where the positivity of the measures is not required. So we 
will be dealing with signed measures. 

Definition 3.1. The modified equilibrium problem asks for two signed measures [i\ 
and [ii minimizing the energy (|1.6[) among all signed measures with supp(/ii) C 
[a, 0], supp(/i 2 ) C [0, 1] and J dpi = |, J" dfi 2 = \- 

The modified equilibrium measures \x\ and \i 2 are unique, and have full sup- 
ports [a,0] and [0,1], respectively. In fact, the modification of the equilibrium 
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Figure 5. A sketch of the modified equilibrium densities in the 
three cases a < — 1, a = 1 and a > — 1. The modified equilibrium 
density becomes negative near on the larger of the two intervals. 

problem comes down to forcing the equilibrium measures to have full supports, at 
the expense of losing positivity of one of the measures near 0. 

In the symmetric case a = — 1 the modified equilibrium measures coincide 
with the usual equilibrium measures. In this case the equilibrium densities ipi and 
ip2 are positive on the whole supports, and around they blow up like an inverse 
cube root. 

In the general case a ^ -1 the equilibrium density on the smaller interval is 
positive on the full interval. The density on the larger interval becomes negative 
in an interval between and xq where xq depends on a in such a way that 

x = so (a) = ^^g^ + 0{a + l) 4 as a -> -1, (3.1) 

see (|3.11l) below. For o^-l both densities il>i(x) and if>2{x) behave like x~ 2 / 3 as 
x — > 0. A sketch of the densities is given in Figure [5] 

Define the logarithmic potentials of a (signed) measure n by 

U"(x) = f log dfi(y), xeC. (3.2) 

J \x-y\ 

The Euler-Lagrange variational conditions [141 133) for [i\ and /12 then say that 
there exist constants l\, I2 G M such that 

2U^(x) + U^(x) =h, for x G [a, 01, 

(3 31 

U^{x)+2U^{x) = l 2 , for ie [0,1], v ' ' 

and these conditions characterize the modified equilibrium measures. For the non- 
modified equilibrium measures we would have an inequality instead of equality 
for x in the gap of the supports. The fact that for \x\ and /12 the Euler-Lagrange 
variational conditions have such a simple form on the full intervals [a, 0] and [0, 1] 
will be important for the further analysis. 

3.2. Riemann surface 

The modified equilibrium problem is easiest to analyze by means of an appropriate 
three-sheeted Riemann surface 1Z. Define 1Z by taking three copies of the Riemann 
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Figure 6. The sheets of the Riemann surface 1Z 



sphere C with cuts 

n :=C\[a,l], fti:=C\[a,0], ft 2 := C \ [0, 1], (3.4) 

and gluing them together along these cuts in the usual crosswise manner, see 
Figure [6] 

The Riemann surface has genus zero and can be defined by the polynomial 
equation 

4a£ 3 - 2(a + l)z£ 3 - 3(a - l)z£ 2 + (a - l)z = 0. (3.5) 
Solving for z, we find a rational function 

Z = Z(0 = 2(a + l)? + 3(!-l)e-(a-iy (3 ' 6) 

which defines a conformal map from £ 6 C to z € 7Z, so that the branch points 
a, 0, and 1 of 1Z correspond to £ = —1,0,1, respectively. The restriction of the 
inverse mapping of (|3 . 6[) to the sheet TZi is denoted by £j 



(i-.Ui^ C. 

The ^-functions map the sheets of 1Z to certain domains 

of C. Then 1Z\ and TZ2 are bounded, while TZq is unbounded. We use 7^ and 7^ 
to denote the arcs bounding TZ\ and 72.2 with clockwise orientation as in Figure [8] 
below. 

3.3. Properties of the modified equilibrium problem 

We can make use of the Riemann surface to prove the following properties of 
the modified equilibrium measures. We will not give all details in the following 
calculations. 
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The modified equilibrium measures are described in terms of the algebraic 
equation 



C Az{z-a){z-l) C 4z 2 (z-a)(z- 1) ° (3 ' ?) 
where z* = z*(a) is a certain solution of 

64(z*) 3 - 48(a + l)(z*) 2 - (15a 2 - 78a + 15) z* - (a + l) 3 = 0. (3.8) 

It can be shown that (|3.8p has three distinct real solutions if a < 0. We use z*(a) 
to denote the middle one of the three solutions and this is the value that is used 
in (33) . 

Proposition 3.2. For a < 0, the following hold. 
(a) The three solutions of (13.71) are given by 



ClW = -/°*M ±t , (3.9) 

J a Z ~X 

{ 2 (z) = -[ 1 ^d X , 
Jq z x 

where ipi ond "02 o^e the densities of the modified equilibrium measures. 

(b) The densities satisfy 

01 (x) = ^— (&,+ {x) - Ci -0*0) , x e (o,0), 

2 ™ (3.10) 
V> 2 (z) = — (Ca,+(x) - C 2 ,-(z)) , i e (0, 1). 

(c) There is Xq = xq(o) G (a, 1) iwira f/ie same sign as a + 1 swc/i £/ia£ 

• if — 1 < a < 0, i/ien ipi( x ) > / or ^ S ( a J 0) arlc ^ ip2(%) < if and only 
if0<x<xo, 

• if a < —1, then il)%(x) > /or x G (0, 1) and ipi{ x ) < if and onZy i/ 
xo < a; < 0. 

See a/so Figure\5\ 

(d) We aawe 

fa + IV 3 

so(o) = 1AQ +0(a + l) 4 <wa->-l. (3.11) 
1U8 

Proof. Let Q, j — 0, 1, 2 be defined by (|3.9j) . so that we clearly have 

Co(*) + Ci(*) + C2(«)=0 (3.12) 

It follows from the variational conditions (|3.3|) that £o,+ = Ci,- on ( a >0) and 
Co,- = C 2> - on (0, 1). 
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Thus if we consider Q as a function denned on the sheet TZj for j = 0, 1, 2, 
then this function extends to a meromorphic function on 1Z. Since (due to the 
normalization J d/ii = J d/12 = 1/2) 

Co(z) = ^ 1 +0(^ 2 ), Q{z) = ~z- 1 + 0(z- 2 ), (3.13) 

as 2 — > oo, the meromorphic function has simple zeros at the three points at 
infinity. There are simple poles at a and 1 and a possible double pole at 0. In 
addition there is a fourth simple zero at a point z*. 

Then the product C0C1C2 is a rational function in the complex plane with a 
zero at z*, simple poles at a, — 1, a double pole at 0, and it behaves as \z~ z as 
z — > 00. This means that 

= 4z2{z Z S a ;l_ iy (3-14) 

Similar considerations show that 

Co(*)&(*) + + &(*X2(*) - - teU-a)U-l) (3A5) 
for some q. Thus Q, j = 1, 2, 3 are the three solutions of the algebraic equation 
f3 3z - q z-z* 

^ 4z(z-a)(z- iy 4z 2 {z~a){z- 1) 1 J 

Inserting ( = Ci(z) = - 1 /( 2z ) + 0{ z ~ 2 ) into (|3T6| shows that q = 2z* + a + 1, 
which gives us the equation (|3.7j) . 

The discriminant of (|3.7I) with respect to C has the form 

Qa(«) 



16z 4 (z-a) 3 (z-l) 3 

where Qi(z) is a certain quadratic polynomial in z that we calculated with Maple. 
The poles a, and 1 of the discriminant correspond to the branch points of the 
Riemann surface. The quadratic polynomial should have a double zero, since oth- 
erwise there would be more branch points. This leads to a condition on z* 1 which 
turns out to be given by (I3.8[) . Again we made these calculations with Maple. This 
proves part (a) of the proposition. 

The relevant solution z* of (|3.8j) is the one that is for a = — 1. This solution 
is then well-defined as a real analytic function for a £ (— 00, 0). We have 

^(«) = -^^ + 0((« + l) 4 ) asa^-1, (3.17) 

which can be obtained from (|3.8I) . The double root of Q%{z) turns out to be equal 
to 

(1 + a) 3 + (6a 2 - 42a + 6)z*(a) - 15(1 + a)(z*(a)) 2 + 8(z*(a)) 3 



x (a) 



18 (1 - a + a 2 - 2(1 + a)z*{a) + (z*(a)) 2 ) 
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which can be shown to also satisfy a cubic equation 

(27a 2 - 46a + 27)a$ - 3(a + l)(9a 2 - 14a + 9)xq 

+ 3a(lla 2 - 14a + ll)x - a(a + l) 3 = 0. (3.18) 

There are three real distinct solutions of (|3.18j) if a < and xo(a) is the middle 
one. The expansion (|3.11[) follows from (|3.18j) and part (d) follows. 

Part (b) follows immediately from part (a) and the Sokhotskii-Plemelj formu- 
las that tell us how to recover the density of a measure from its Cauchy transform. 

Finally, to prove part (c), we suppose that — 1 < a < 0. It can then be shown 
from the above formulas (it is not immediate, however) that < xq(o) < 1. Since 
xo(a) is a zero of the discriminant, the cubic equation (|3.7j) has a double solution 
if z = Xo(a). Since < xo(a) < 1, we have £o t +(xo(a)) = £,2,+ (xo(a)) and £i(xo(a)) 
is real. Thus we have ^o,+ (xo(a)) = ^2,+ (xo(a)), which means since £o,+ = £2,- 
on (0, 1), that V>2 vanishes at Xo(a) by (|3.10p . Since Xo(a) is the only zero of the 
discriminant in (a, 1), it also follows that Xo(a) is the only zero of ip2, and that tpi 
has no zeros. Thus ipi > on (a, 0). It is a consequence of the fact that the point 
of the Riemann surface is a double pole, that ipi and ip2 have opposite signs near 
0. Thus ip2(x) < for < x < xo(a) and part (c) of the proposition is proved in 
case — 1 < a < 0. 

The proof for a < — 1 is similar. □ 



4. Steepest descent analysis of the RH problem 
4.1. First transformation 

We start from the RH problem for Y with m = nz = n. We also take a < 
close to —1 but for the moment it is arbitrary and fixed. We use the modified 
equilibrium measures \i\ and /12 that are supported on the two intervals [a, 0] and 
[0,1] respectively. 

Define ^-functions by 



where we use the main branch of the logarithm. Hence g\ is defined with a branch 
cut on (—00, 0] and #2 with a branch cut on (— 00, 1]. The boundary values of the 
g-functions along the real axis are given by 




(4.1) 




(4.2) 
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were U^ 1 and U^ 2 are the logarithmic potentials ()3.2|) . From (|3.3j) and (|4.2[) we 
obtain 

9i,+ (x) + 9i,-{x) + g 2 ,±(x) = -h± -ni, x € [a,Q], 
9i(x) +g 2 ,+(x) +g 2 -(x) = ~h, x G [0,1], 

Now define the first transformation 7 H- J 1 as 

/l 

T(z)= o e - 2 "(^ + 5-) | 
\0 e~ 2 ™' 2 

/ e -2n(g 1 (z)+g 2 (z)) q q 

X g2n(gi(2)+Zi+27ri) g I . (4.4) 

^ e 2n(s 2 (z)+«2) 

This transformation normalizes the RH problem at oo, since the g- functions 
behave like \ logz + 0{z~ x ) as z -> oo. Thus T(z) = 7 + 0(z- 1 ) as z — > oo. The 
jumps for T are conveniently expressed in terms of the two functions tpj, j = 1,2 
defined by 



ipi(z) := -2gi(z) - g 2 (z) ~h + 



Vi{z) ■= -252(2) - 91(2) - Z 2 + 



^7ri for Imz > 0, 
— for Im z < 0, 

iri for Im z > 0, 

— 7ri for Im z < 0. 



(4.5) 



Then by (|2.1[) , (|4.4p . and (|4. 3[) one obtains the jump matrices 

/ e 2n Vl . + (x) w ^ q\ 

J T (x) = e 2n<pi,-(*) for a; € (o,0), (4.6) 

V 1/ 

J T (ar)= 10 for x G (0,1). (4.7) 

V e 2n ^-^J 

Thus T satishes the following RH problem. 

• T is an analytic 3x3 matrix valued function on C \ [a, 1] , 

• T satisfies jump conditions T + (x) = T^(x)Jt(x) for x G (a, 0) U (0, 1), with 
J T given by ([476]) and ([4~7]) . 

• T is normalized at infinity: 

T(z) = 1 + (z~ x ) as z ^ oo, (4.8) 

• near a, and 1 the function T has the same behavior as Y, see (|2.3p . (|2.4I) . 
and ([2751. 
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Figure 7. The contour Eg and the lenses around [a, 0] and [0, 1]. 



4.2. Second transformation: opening of the lenses 

The functions e 2nVk - ± 1 k = 1,2 appearing in the jump matrices Jt in (|4.6|) and 
(|4.7[) are rapidly oscillating for large n, since the boundary values <pk,± are purely 
imaginary. One may easily check that 



ipi ±(x) — ^f2ni / dfXi(s) for #6 (a, 0), 



<P2,±{%) = ±27ri / dfi 2 (s) for x € (0, 1). 
Jo 



(4.9) 
(4.10) 



The oscillations are turned into exponential decay by the so-called opening 
of the lenses. Choose smooth paths £f connecting a and with E+ in the upper 
half-plane and Ej~ in the lower half plane. Similarly, choose paths Ej" and 
connecting and 1. Define 



E s := Sf U U [a,l]. 



(4.11) 



The intervals [a, 0], [0,1] and the paths T,f define 4 bounded regions that arc 
referred to as the lenses around [a, 0] and [0, 1], see Figure [7] 

For j = 1, 2, let Vj be a simply connected neighborhood of Aj such that the 
analytic factor hj in the weight function Wj is analytic and non-zero in Vj. We 
assume that E^ C for j = 1,2. Then iUj has an analytic continuation from Aj 
to Vj with some cuts, which we also denote by wj\ 



Wl (z) = {z-a) a {-zfhx{z), for z € Vi \ ((-co, a] U [0,oo)) . 
w 2 (z) = /(l - z) 7 /i 2 (^), for z e V2 \ ((-oo, 0] U [1, oo)) . 



(4.12) 
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Then, following [29j section 4], we define the next transformation T i— > S by 
/ 1 0^ 



S(z) 



S(z) 



T(z) 



T(z) 



T{z) 



T{z) 







1, 
1 

w 1 {z)- 1 e 2n ^ zS > 1 
1, 



1 N 

1 

\-w 2 {z)- l e 2n ^ z ) 1, 
1 N 

1 



K w 2 (z)~ 



-Vf2(z) 1 



z in the upper part of 
the lens around [a, 0], 

z in the lower part of 
the lens around [a, 0], 

z in the upper part of 
the lens around [0,1], 

z in the lower part of 
the lens around [0, 1], 



(4.13) 



(4.14) 



and 



(4.16) 



S(z) = T{z) elsewhere. (4.15) 

It is clear that this transformation does not affect the behavior at infinity. 
The jump matrix J$ for S on the intervals (a, 0) and (0, 1) are 

' / wi{x) 0\ 

-w^x)" 1 , x € (o,0), 
= < \ 1/ 

/ w 2 {x)\ 

1 , x G (0,1). 

[ V-waCa;)- 1 / 

The transformation has introduced jumps on which are 

1 o\ 

1/ 

1 0\ 
1 , z E T,f. 

\w 2 {z)~ 1 e 2n ' f ' 2 ^ 1/ 

Finally, the behavior near a, and 1 changes because of the factors wj 1 in 
the transformation (|4.13[) - (|4.14[) . The Riemann-Hilbert problem for S then reads: 

• S is analytic on C \ Sg, 

• S 1 has jumps S + — S-Js, where Js is given by (|4. 16[) and (|4.17[) . 

• S(z) = I + (z^ 1 ) as z -> oo, 



(4.17) 
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• near the endpoints of the intervals S behaves as 
S(z) = O 



e x (z) e 2 (z) 
ei(z) e 2 {z) 
ei(z) e 2 (z) 



(4.18) 



where 



ei(z),e 2 (z) 




S(z) = O 



where 



ei(z),e 2 (z) 




f ei(z) e 2 (z) 
S{z) = 0\e 1 {z) e 2 (z) 
,ei(z) e 2 (z) 



if a < 0, 
if a = 0, 

if a > 0, z inside the lens, 
if a > 0, z outside the lens, 

as z — > 1, 



if 7 < 0, 
if 7 = 0, 

if 7 > 0, z inside the lens, 
if 7 > 0, z outside the lens, 



as z 



(4.19) 




(4.20) 



where 



ei(z),e 2 (z) 




if j8 < 0, 
if j3 = 0, 

if /3 > 0, z inside the lenses, 
if /3 > 0, z outside the lenses. 



For later analysis it will be important to know how Re (pi and Re (p 2 behave 
on the lips of the lenses. From (I4.17|) we see that we would like to have 



Re tpj(z) < for z e T,-- 



(4.21) 



for j = 1,2. The inequality (|4.21j) will indeed hold if /ij is a positive measure, and 
this can be proven using the Cauchy-Riemann equations. Recall however, that /ii 
and \x 2 are signed measures. The inequality (|4.21[) will be violated for z on the 
parts of that are close to the interval where [ij is negative. By Proposition (I3.2|) 
we have that Hi is negative on (xo(a),0) if a < — 1, and that \x 2 is negative near 
(0, x (a)) if -1 < a < 0, where x (a) = 0((a + l) 3 ) as a -> -1. 
We write (pj(z;a) to emphasize the dependence on a. 
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Lemma 4.1. There exist positive constants Co and C\, independent of a, such that 
for every a sufficiently close to —1, we have 

Reip 3 {z;a) <C a \a+l\\z\^ 3 -dlz^ 3 , z G S^, \z\ < 1/2. (4.22) 

for j = 1,2, 

Proof. We have by the definitions (|3J]) . (T4~T]) . and (j4~5j) that ^ = £? - Co for 
j = 1, 2. The constant in (|4.5p is taken so that ^ (0) which means that 

ipj(z;a)= / (C 3 (s; a) - ( Q (s; a))ds (4.23) 



where Co , Ci ; C2 are the three solution of the cubic equation (|3.7|) , where we em- 
phasize the dependence on a. 

As s — > with Im s > we can compute from (|3.7[) that 



Co(s; a) = c q lu 2 S - 2/3 + Ci los' 1/3 + 0(1), (4.24) 

&(«; a) = cos" 2 / 3 + ds- 1 ' 3 + 0(1), (4.25) 

C 2 (s; a) = cqlus- 2 ' 3 + Ci lo 2 s- 1/3 + 0(1) (4.26) 

uniformly for a close to —1, with real constants Co = (— z*(a)/(Aa)) 1 / 3 and c\ — 
— (2z*(a) + a + l)/(12aco). Because of the behavior (|3.17|) of z*(a) we have that 

co = co(a) = -^-(a + 1) + 0((a + l) 2 ) 

ni/3 (4-27) 
ci =ci(a) = -^- + 0(a + l) 

as a —> — 1. 

Using (j434l) in (|4T23|) we find 

<p 2 {z; a) = 3 3/2 zc (a)z 1/3 - ±3 3/2 i Cl (a)z 2 / 3 + 0(z) (4.28) 

as z — > with Imz > 0. By (|4.27|) we have that ci(a) tends to a negative constant 
as a — ¥ — 1. Since we may assume that the lens is opened with a positive angle at 
0, we find that 

Re ( -IVZic.Mz 2 ' 3 1 < -CM 2 / 3 , : ^ 



;V3 lCl (a)z 2 / 3 ) < -d|z| 2 / 3 



2 



for some constant C\ > independent of a. Using this in (|4.28j) we obtain (|4.23j) 
for j = 2 and z e in a fixed size neighborhood of z — 0, say \z\ < ro- The 
inequality (|4.23|) then also holds for \z\ < 1/2 (maybe with different constant 
Ci), since <p%(zya) — > (p2(z;— 1) as a — > — 1 uniformly for ro < \z\ < 1/2, and 
Re(f2(z;— 1) < -C3 < for z e S^, r < \z\ < 1/2, and some C 3 > 0. 

The inequality (|4.23l) for j — 2 and z£S 2 and for j = 1 follow in a similar 
way. □ 

It follows from (|4.22[) that we indeed have that Re tpj(z; a) < for z € £^ , 
except for z in a small exceptional neighborhood of z = 0, whose radius shrinks as 
((a + l) 3 ) as a -> -1. 



24 



Klaas Deschout and Arno B.J. Kuijlaars 



4.3. Outer parametrix 

The next step is to construct an approximation to S. This so-called parametrix 
consists of an outer parametrix N that gives an approximation away from the 
cndpoints a, and 1 and local parametrices P around each of the endpoints. 

4.3.1. Riemann-Hilbert problem for N. The outer parametrix N should satisfy 
• N is analytic on C \ [a, 1] , 



• N satisfies the jump conditions 

2V_(ar) 



N+(x) = < 





~JOl(l) _1 




N—(x) [ 

-w 2 (x) _1 

N is normalized at infinity: 

N(z) =J+0 



wi(x) 

1 

w 2 (x) 

1 




if a; € (a, 0), 
if x G (0,1), 



(4.29) 



— | as z 



(4.30) 



4.3.2. Solution in a special case. First we will find a solution N to this problem for 
the case that both weights w\ and are identically 1 on their respective intervals. 
We solve the problem for N by using the Riemann surface 1Z introduced before. 
There is a similar construction in [9] and so we do not go into much detail here. 

The function £ = £j(z) maps the sheet TZj of the Riemann surface onto the 
domain TZj as shown in Figure [5] that are separated by two closed contours 71 and 
72 that we orient in the clockwise direction. We write 

Pj =Pj(a) =^(00), j = 0,1,2. 
A solution N is given in the form 

(noMz)) ivo(aw) ivb(6W)\ 

N(z) = NiMz)) NiW)) , (4-31) 



with the following functions iVj that are analytic on C \ (7^ U 7^ 



pjO^-i) 1 / 2 



»=0 



1)1/2 1J 



(4.32) 



with appropriate modifications if po — 00 (which happens if a= — 1). The branch 
cut of the square roots (£ 2 — l) 1 / 2 and (p 2 — l) 1 / 2 in (|4.32p is defined along 7j t ~U7^~, 
see Figure |S] It can then be checked that the matrix function N defined by (|4.31[) 
(|4.32|) indeed satisfies the conditions in the Riemann-Hilbert problem for N. 
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Figure 8. The images of the sheets of 1Z under the functions £j. 
The location of pi — £1(00) and P2 — £2(00) depend on the choice 
of a, but the contours 7^ and 7^ do not. 



4.3.3. Solution in general case. Now we turn to the problem for N, with jumps 
involving W\ and wi. This can be solved using analogues to the Szego function as 
in [37J US]- We look for three functions .Do, D\ and D2 satisfying 

(a) D$ is analytic and non-zero in C \ [a, 1], 

(b) D\ is analytic and non-zero in C \ [a, 0], 

(c) D2 is analytic and non-zero in C \ [0, 1], 

(d) Dq, D\ and D2 have limiting values on (a, 0) and (0, 1) such that 

Di.+ Di_ 



Dn 



and 



D 



2,+ 



Do 



Do 



w\ on (a, 0) 



N(z) 



Do,- A),+ 
Having Do, D\ and D2 we define N by 

L 

^Di(oo)- 1 




W2 on (0, 1). 





We will abbreviate this as 



N(z) = D^N{z)D(z), 



(4.34) 
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and one can check that N indeed solves the Riemann-Hilbcrt problem for N. 
From the jump properties of Dq, D\ and D 2 it follows that 

(A,£»iA>)+ = {DqDxD^)- on (a, 0) and (0, 1). 

Thus D0D1D2 is analytic across these cuts. We also make sure that the possible 
singularities at a, and 1 are removable. Then DqDiD 2 is a constant and we can 
choose a normalization such that 

D Q D X D 2 = 1. 

4.3.4. Szego functions. In order to hnd Dq, D\ and D 2 we write 

D j (z)=V(^ j (z)), i = 0,l,2, (4.35) 

for some yet to be determined function T> on the £-Riemann sphere. Recall that 
£0, £1 and £2 are the mapping functions from the respective sheets of the Riemann 
surface TZ to the Riemann sphere. 
Then V has to satisfy 

(a) T> : C \ (71 U 72) — > C is analytic and non-zero. 

(b) On 71 U 72 there is a jump 

V + (£) = w j (z)V^), £e 7 j, j = l,2, (4.36) 
where z = z{£) is related to £ by (|3.6|) . 
Then by taking logarithms we get 

log© + (0=log«; j («)+logl>_(0, ^7+ 3 = 1,2, (4.37) 
which by the Sokhotskii Plemelj formula is solved by the Cauchy transforms 

log2>(0 = J- ( [ l0gWl{ f }) d S+ [ l0gMz{s)) ds) + C\ 
2m \J 71 s - £ J J2 s - £ J 

where C\ is an arbitrary constant. Thus 
X>(£) = Cexp 

with C = e Cl and 
Dj(z) = Cexp 



1 ( f logt«i(z(s)) d f log w 2 (z(s)) ^ 



1 f /" logwi(z(s)) rfs+ /• log w 2 (z(s)) ^ 



for j = 0, 1, 2. The constant C can be taken so that DqD\D 2 = 1. This completes 
the construction of Do, Z?i and £>2 and therefore of N. 

Example. In the case where hi = 1, /12 = 1, we can evaluate the Szego functions 
explicitly. Indeed, we find for 

Wl (x) = w[ a ^\x) = (x- a) a (-xf, x G (o,0), 

« , (4.39) 
W2 (x) = w f n \x) = a^(l - xf, x g (0, 1), 
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that 



D (z) = C 



fr(*)-6(c») \°Y 6>(*)-6M 

(eo(z)-Ci(oo))(£o(*)-&(oo ))V 



^ ) = C r fl) 6W-6(a)J U(*)-6(i)J 
v (, ^i(*)-fr(oo))(fr(*)-fr(°°)) V 



(4.40) 



e 2 (z)-6(a)y V v 

(6H-6M)(6W-6M)\^ 



with appropriate choice of branches for the exponents. 

4.3.5. Behavior of Szego functions near 0. From (|3.5j) and the choice of branches 
£o, fi, £2, we obtain 

fo(«) = -C(o)w^V3 + ( i! ) i 

£1(2) = -C(a)z 1/3 + C(z), for±lmz>0 (4.41) 

£2(2) = -CCa^^ + O^), 



with 



1 \ 1/3 
a — 1 x 



We use this in (|4.40jl together with ^(a) = -1, £2(1) = 1, £1(00) = pi(a), 
£2(00) = p 2 (a), to obtain the leading behavior of the Szego functions at for the 
case hi = 1, /i 2 = 1. It follows from (I4.40j) that 

D (z) = C (- Pl (a)r P2 (ar ( ~ P ffiffi (o) )" ^^(l + 0{z^)), 
D 1{Z ) ^ C (a Pl (a)T P2 {ay ( "^^"^ ^(l + O^ 3 )), 
D 2 (^) = C(-^(a))>2(a)^ ( ~^) 2 2(a) ) ' e*™/V"(l + OfrV)) 

(4.42) 

as z — > with ±Imz > 0. 

The effect of the analytic factors h\ and /i 2 comes in the form of contour 
integrals 

1 f hghj(z(s)) 



2%ih '-ds, 3 = 1,2, ^C\ 7 ,, 
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see (|4.38j) . Because of analyticity we can deform jj to a contour t| in the region 
IZj , which leaves the integral unchanged if £ € C \ IZj and picks up a residue 
contribution of \og(hj(z(£)) in case £ € 7£j is close to jj, in particular if £ is close 
to 0. In this way we find the following behavior as £ — > 0, 

1 / log^(z(,s)) rfs= |c i + O(0, £GC\7^ (443) 



27ri/ 7J s-£ [ Cj -log^(0) + O(0), £eft 3 

where 

1 / log/i j (z(s))_ dg 



By the change of variables z(s) = x, s = £j (z), we turn this integral into an integral 
on a counter that circles around Aj in counterclockwise direction. Bringing this 
integral to Aj we obtain 

^ = ^/ 4j 1<>6ft ' (l) ((l)-"(§) + ) WdI ' (4M) 

Combining (|4.42l) . (|4.43l) . (|4.44|) we find that for general analytic factors we 

have 

D (z) = Ce<*+<* (-pi(a)r P2 (ay ( "^^ l^^V^Cl + <^ 1/3 )), 
Dx(z) = Cfe 1 (0)e^+ e » (apxfo))^^ ( z^^ )* ^/3 (1 + O(z^)), 

D 2 (z) = Ch 2 (0)e c ^ {- Pl {a)TP2{aV f 1 ^^) /3 (1 + 0(z 1/3 )) 

(4.45) 

as z — > with ilmz > 0, with c% and C2 given by (|4.44|) . 

For a = — 1 we have — pi(a) = p 2 (a) = 775 arid C(a) — 2~ 1 / 3 . Since all 
quantities depend analytically on a, we find from (|4.45l) 

D (z; a) = Ce c ^ - , J|^ e ±W 2 -2/3/3 (1 + 0(z i/ 3) + c(a + 

D x (z; a) = C/ ll (0)e Cl+C2 -Jg^ ^l + 0(^/3) + 0{a + ( 4 .46) 

2I/ 3 

D 2 (z; a) = Cfe 2 (0)e Cl+C2 n - - e T/3 " /3 z^ /3 (l + ©(z 1 / 3 ) + 0(o + 1)) 
as z — » and a — > — 1. 

4.3.6. Behavior of N around the endpoints. For the further analysis of the Riemann- 
Hilbert problem for S we need to know the behavior of N around the endpoints 
a, and 1. To that end we also need to know how the Szego- functions Do, D\ and 
D 2 behave around these points. 
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The functions Z?o, 1,2(2; Wi, w 2 ) are multiplicative in w\ and w 2 . Then we can 
split off the analytical factors hi, h% from the weights w±, w 2 and write 

Dj(z;wi,w 2 ) = D j (z;w[ a ' l3 \w 2 '^)D j (z;hi,h 2 ) 



) associated with the weights (|4.39[) are 



where the functions Dj(z;wi"' ,wi^ 
given in (14.401) above. 

The explicit expressions allow us to prove the following proposition: 



Proposition 4.2. Around the branch points N has the following behavior: 

A 1 + 2q . , 2a-l \ 

z — a) * (z — a) * 1 \ 

O (z — a) {z — a)^r~ 1 as z — >• a, 

\(z — a)~^^~ (z — a)^r~ lj 



N(z) = < 



O 



O 



1+27 
4 

1 + 2^ 



f(z - I)" 

(z-1)- 

\(z-l) 4 

/ 20 + 1 g-1 
/ Z 3 2 3 



1 (z-iy 

1 (z-ir 
1 (z-l) ; 



£-1 . 

Z 3 

2<3 + l g-1 f3-l 

Z 3 2 3 Z 3 

23 + 1 g-1 
Z 3 2 3 Z 3 



1. 



as z — > 0. 



(4.47) 



Proof. From the expressions (|4.31[) and (|4.32[) together with the behavior of the 
mapping functions £o,£i,£2 around the branch points we find ()4.47|) for the case 
wi = 1, W2 = 1, (in which case of course a = /3 = 7 = 0). 

For the general case we first note that the functions DAz\hi,h-2) remain 
bounded and bounded away from for analytic and non-zero hi and h 2 . For the 
Szego functions (|4.40[) associated with the pure Jacobi weights we have 



D (z;w[^\4^) Di(z;w{^\ W i^) D 2 (z; W {^\ w ^) 
' O ( \z - a)- a ' 2 (z~a) a / 2 l) 



( (z - 1)-t/ 2 1 (z-l)^/ 2 j 

Q ( z -2f3/3 z 0/3 z P/3^j 



as z — > a, 
as z — > 1, 



as z 



0. 



where we use that 



£j(z) = C(z 1/3 ) asz^O, 
ij{z) = 0((z- a) 1 ' 2 ) as z - 
i ] {z) = 0{{z-l) 1 ' 2 ) asz- 



for J =0,1,2, 
a, for j = 0, 1, 
1, for j =0,2, 



while £i(z) is analytic around z = 1 and (, 2 {z) is analytic around z 



(4.48) 



□ 
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4.3.7. Symmetries in the outer parametrix. In this subsection we give two sym- 
metries in the functions N that will be useful later on. Recall that ./V is the outer 
parametrix in the case that wi and W2 are identically one, see subsection ([4.3.21) . 
A first symmetry deals with the inverse of N. 



Proposition 4.3. For every a < we have 

N-\z) = N T (z), zeC\[a,l]. 

Proof. Define X by 

X(z) := N(z)N T (z) zeC\[a,l] 



Using the fact = ( J 



'Nl 



we find that on (a, 0) and (0, 1), 



xz 1 x+ = nz t nz 1 n + nJ = nz t j r nI 



N~ 



J 



NJ 



Nl 



N 

+ N 



(4.49) 
(4.50) 

(4.51) 



Since N(z) tends to I as z — > oo we have X(z) = I + (i) as z — > oo. 

By the behavior of N near the branch points (see (|4.47l) for the case a = /3 = 
7 = 0) we obtain that X has no poles in a, or 1, and we conclude by Liouville's 
theorem that X(z) = I everywhere and (|4.49[) follows. □ 

As a corollary of this proposition and the expression for the Dj in (|4.38|) we 
then also find the behavior of TV -1 around the branch points, since from (|4.34j) 
and (j4^9l) 



N- 1 (z) = D(z)- 1 N(z) 1 D c 



(4.52) 



We then obtain 



N~ 




(z — a)' 
[z-a)~ 
1 

(z-lf 
1 

(z-ir 



as z — > a, 



as z — > 1, 



3 


z 


3 


z 


3 


as z 


/3 + 1 
3 


z 


P+l 
3 


z 


3 1 





0. 



(4.53) 

A second symmetry relates the functions N for different values of a to each 
other. We use N(-;a) to denote the dependence on a < 0. Let a be the Mobius 
transformation leaving and 1 invariant and mapping —1 to a, i.e., 

2az , , , (1 — a)z 



(a + l)z + a — 1 



(a + l)z - 2a 



(4.54) 
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Proposition 4.4. For a < 0, a ^ — 1 and any z E C \ [a, 1] we have 



N(z;a)=N(j-^;-lj N {6-\z); -l) (4.55) 

Proof. The proof is similar to the proof of the previous proposition. We do not 

give details. □ 

4.3.8. Behavior at z = 0. In the next section we need the leading term for N(z; a) 
as z — > 0. 

Lemma 4.5. As z — > we /icrae 

AT(z;a) = iV (a)z^ +0(1) (4.56) 



—uj 1 w 2 I /or Im z > 0, 

2 i V I (4 ' 57) 

ur 1 w /or Imz < 0. 



Proof. By Proposition 14.41 and the fact that 8 a 1 (z) = z + 0(z 2 ) it suffices to 
compute the leading term of N(z; — 1) as z — > 0. We use the factorization (where 
we suppress the argument z on the right hand side) 



iV(z;-l) = diagfl EllA^H ntf-D 1 ") 

\ pi— P2 P2—P1 J 



'(& -pi) (Co -vi) (6 -pi)(6 -P2) (6 -pi)(6 -P2) N 

Co ~ P2 £l - P2 6 - P2 

Co - pi £i - pi 6 - Pi 

A careful analysis of all the functions and constants involved then shows that for 
Im z > 



32V6 

and for Imz < 

1 



iV(z;-l) = — I 1 I (-w 1 a; 2 ) +0(1) as z ^ 0, (4.59) 



7V(z;-l) = ^ T7¥ z"3 | i | ( w a ! w ) + G (l) as z ^ 0. (4.60) 



Together with (|4.55l) this proves the lemma. □ 
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4.4. Local parametrices 

4.4.1. Local parametrices around ±1. The outer parametrix TV is intended as an 
approximation to S. However the approximation cannot be good around the branch 
points a, and 1. Indeed, the entries of S(z)N~ 1 (z) will typically diverge as z tends 
to one of the branch points. 

The solution to this problem is building local approximations around the 
branch points, called local parametrices. The appropriate construction around a 
and 1 is standard, and uses the Bessel model parametrix as defined in [29], equa- 
tions (6.23)-(6.25). Let U-i and U\ be disks around respectively —1 and 1 of fixed 
but small enough radius: U\ should be contained in V%, such that u>2 is well-defined 
on Ui \ [l,+oo). Similarly we must have f7_i C Vi, and additionally a must lie 
inside U-i. On these disks we construct 3x3 matrix valued functions P_i and 
Pi that satisfy the same jumps as S, see (|4.16|) - (|4.17|) . and match with N on the 
boundary of the disks: 

P-i{z)N(z)^ 1 = I + O [ — ) for z G dU-i, 

(4.61) 

for z £ dUi, 



P l {z)N{z)- 1 =I + o{^j 



as n — » oo. The O-terms are uniform in z. For details of the construction of P_i 
and Pi we refer to [29], where the Bessel model parametrix was introduced, and 
[31j . where it was also used in a 3 x 3 matrix valued Riemann-Hilbert problem. 

4.4.2. Local parametrix around 0: statement. Around we need a new kind of lo- 
cal parametrix Pq on a disk Uq around the origin. There are a number of difficulties 
to obtain the desired matching condition 

P {z)N{z)- 1 = 1 + 0(n- K ) for z £ dU 

with some k > 0, that in fact we are unable to resolve. The best we can do is to 
construct Pq such that Pq(z)N{z)~ 1 remains bounded as n — > oo for zona circle 
of radius that decays like n -1 / 2 asm oo. 

Thus the disk Uq should be shrinking as n increases, and for definiteness we 

take 

radius U := rT^ , (4.62) 

and we assume n to be large enough so that Uo is contained in V\ n Vi- Then 
consider the following Riemann-Hilbert problem for the local parametrix Pq around 
0. 

• Pq is analytic on Uq \ Eg, 

• Pq has jumps 

Po,+ = Po,-Jp onEsflUo, where J P = J s , (4.63) 

see Km - KTT^ . 

• Pq(z) behaves in the same way as S(z) as z — > 0, see (|4.20p . 
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• PoN- 1 



remains bounded on the boundary of Uq, 



P (z)N(z) 



-l 



0(1) for z e dU t 



(4.64) 



as n — > oo, where a = a n depends on n as in (jl.9|) . 

Notice that the matching between Po{z) and N(z) does not improve with 
increasing n. Indeed, the matrix Pq(z)N(z)~ 1 does not tend to / as n — > oo for 
z 6 OUq. The matching (14.64[) is the best we can obtain without modifying the 
outer parametrix N. However, with a = a n as in (|1.9|) we will be able to find a 
3x3 matrix valued function Z n (-; a) such that 



The explicit expression and special properties of Z n will allow us to create, in the 
final transformation, a jump on OUq that tends to / as n — > oo. 

Since the dependence on a will be important, we emphasize that most notions 
depend on a and have limiting values as a — > —1. As before, we will not always 
explicitly indicate the dependence on a, but sometimes we do. 



4.4.3. Reduction to constant jumps. We factor out the tfi and u^-functions from 
the jump matrices (|4.16[) - (|4.17[) . Define for z£[/ the matrix valued functions 




o 



(4.65) 






if 2 (z) - 2(pi(z) 





(4.66) 




(4.67) 



We look for P in the form 



P Q {z) = P {z)e nK ^W- 1 {z). 



(4.68) 
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In order that Pq has the jumps Jp, we should have Po,+ = Po,-Jp with 

1 o\ 

e ±f37ri 1 on T,f n U , 

1/ 

/l o o\ 

1 on EfnU , 



Jp = ( 



( e Pm 0\ 

-e^ on [a,0] D U , 
{ 1/ 

l\ 

1 on [0,1] n U Q . 

-1 0/ 



(4.69) 



4.4.4. Functions f(z) and r(z). Note that the jumps (I4.69|) are exactly the same 
as the ones for see Figure |31 except that the jumps for \E' are on unbounded rays. 
Recall that 9(z;t) also depends on r which appears in the asymptotic condition 

Our aim is to construct Pq of the form 

P (z) =E n (z)v(nif(z);nh(z) 



(4.70) 

where f(z) is a conformal map and t(z) is analytic in Uq. The matrix valued func- 
tion E n (z) is an analytic prefactor, which will be defined in the next subsection. 
We are going to choose f(z) and t(z) such that 



e 



{n^ f(z);n^T(z) S j + nA(z) = for z e U \ ' 



(4.71) 



where & is given by (|2.9|) and A is given by (|4.66j) . When this condition is satisfied 
there will be no exponential growth (as n — > oo) in Po(^) se e (|4.68p . and so there 
is a chance that we can match it with N. 

Define functions Ai(z) and \2(z) on C \ 1 by 



Ai(z) := 



-z 3 ((pi(z) 
-z~r, ((px(z) ■ 



- uj 2 ip 2 (z)) 
UJLp 2 (zj) 



for Imz > 0, 
for Im z < 0, 



(4.72) 



\ 2 {z) 



(4.73) 



— z 3 (<Pi{z) + LJ(f2(z)) for Imz > 0, 
— (<pi(z) + u> 2 ip2{z)) for Imz < 0, 

where ipi and ip 2 are given by (|4.5[) . It can be checked that the functions Ai(z) and 
A2(z) have no jumps on (a, 0) or (0, 1). Since ipi and ip 2 are bounded, Ai and A2 
have analytic continuations to C\ ((—00, a] U [1, 00)). Since the (^-functions depend 
on a, so do the A-functions, and we write A(z; a) to emphasize this fact. 
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It may be checked that Xi(z;a) and \2(z;a) converge uniformly in a neigh- 
borhood of zero as a — > — 1, and also that 

3 • 2 2 / 3 

Ai(0; a) = — (a + 1) + 0(a + l) 2 , 

, (4-74) 

A 2 (0;a) = + 0(a + l) 

as a — > — 1. Then Re A 2 (.z; &) > for a close enough to —1 and z *E Uq, and we can 
define the following analytic functions in a neighborhood of 0. 

Definition 4.6. For a close enough to —1 and z 6 Uq we define 

8 3 

/0) = /0;a) := ^^A 2 (z;a)2, 

r(z) = t(z; a) := -. 

A 2 (z;a)2 

Then / is a conformal map with /(0) = 0, and f(z) is real for real arguments 
z. By (|4~V4|) and (|4~75]) we have as a -> -1. 

/'(0;o) = V^ + 0(o+l), 

1 , (4-76) 

r(0;a) = — (a + 1) + 0(a + l) 2 . 
v2 

Without loss of generality we can now assume that the lips of the lenses are 

chosen such that / maps EgnE/o into X*. Then ^ (n^f(z); n^T(z)^j is well-defined 

and analytic in Uo \ Eg. It remains to check the condition (|4.71l) . By (|2 . 10[) and 
(|4~75) we find: 

9 k {nlf(z);Jr(z)) = _|V/(z)§ - nr(z) W 2fc /(^ 

= -y (wMA2(z)+w 2 MAi(;8)) (4.77) 

By (|4.72[) and (|4.73[) the right hand sides are exactly minus n times the 
components of A, and (|4.71l) follows. 

4.4.5. Prefactor E n (z). Next we define a suitable analytic prefactor E n (z) = 
E n (z;a) such that the local parametrix 

P (z) = E n {z)V (n^f{z);nh(zfj e n ^W{z)- 1 (4.78) 

satisfies the matching condition (|4.64|) with N on 8Uq. Also wc obtain an expression 
for the function Z n in (14.651) . 

Denote by A(z; r) the right-hand side of (12.7[) without the exponential factor 
e c>( z ; T ) anc [ t ne error factor I + 0(z~*). Thus 




A(z;r) := J— e T 'V 1 I Cl±B±, ±lmz>0. (4.79) 
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Then we define 

Definition 4.7. For a close enough to —1 and z € Uq we define 

E n (z; a) = N{z; a)W(z; a)^ 1 (n*f(z; a);n?r{z; a)) . (4.80) 
where f(z\ a) and r(z; a) are given by (|4.75l) . 

Proposition 4.8. The function E n defined by (|4.80|) is analytic on Uq. 

Proof. The function A from (I4.79|) has jumps on the real line, given by A + = A- J a 
with 

/ e^ 1 0\ / 1\ 

J A =\-e^ 1 onR", J A = 10 onR+. (4.81) 
\ 1/ \-l 0/ 

Then A (ni f{z); n^r(z)^ has the corresponding jumps on R n C/q- A straightfor- 
ward calculation reveals that N(z)W(z) has exactly the same jumps on Uq n (a, 0) 
and J7o H (0, 1), and hence E n (z) has no branch cuts in ?7o- 

Also £"„ can have no pole in 0. By (|4.47p and the definition of W (|4.07p we 
have that 

N{z)W{z) = O U^) as z -> (4.82) 
From (|4.79|) and the fact that / is a conformal map with /(0) = we obtain 



.4 



(n*/(«);n*T(*)) * = O (z - ^) (4-83) 



Then by (|4T5D|) and flUgg) and gggl we see that £ n (z) = and therefore 

the isolated singularity at is removable. □ 

We remark that by a similar argument we have that E~ l {z]a) is analytic in 
Uq and in particular at z = as well. 

4.4.6. Matching condition. We show that the matching condition (14.041) holds, and 
we compute Z n from (|4.05[) . 

Proposition 4.9. The parametrix Pq defined by (|4.08p . (I4.70[) . (|4.80[) satisfies the 
matching condition (14.64)) as n — > oo with a = a n as in (|1.9[) . The matching (|4.05[) 
holds with 

1 T(z;a)(nT 2 {z;a) + 9(3) 
A n{z, a) — , 3 , 

oi ■ i '° Z3j(z;a)3 

x D- 1 I 1 I (V2« 1 -l)Doo. (4.84) 
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Proof. The local parametrix Pq, written in full, is given by 

P (z) = N(z)W{z)A (n% f{z);n%T{z)\ ~ * (n^f(z), nK(zj\ e nK{z] W{z)-\ 

(4.85) 

where all functions also depend on a. 

We first observe that n^r(z; a) remains bounded for z € Uq as n — > oo. Here 
we need the fact that we took the radius of U$ to be n _ 2 7 and that a = a n = 
-1 + C(n _1 / 2 ). Since t(0; -1) = 0, see (|4J6]l we indeed obtain 

T (z;a) = t(z;-1) + O (a - a„) = 0(n- 1/2 ) for \z\ = n - 1/2 , (4.86) 

as n — > oo. 

It follows that we can use the asymptotic expansion (|2.22j) for ^(z;t) in 
(|4.85p . since (|2.22p is uniformly valid for r- values in a bounded set. Then by 
combining (|2.22p and the definition (|4.79p of A we get that a number of factors 
cancel, and what remains from (|4.85[) is 

Poiz)N(z)- 1 = N(z)W(z)B ± 1 

x ( /+ ( * l} f ( ;^ )} +0{n-i)) B±W{z)-iN(z)-* (4.87) 

as n — > oo uniformly for z € 8Uq. 

For z € dUo we have \z\ = and the entries of NW are 0{z^~3~) asz->0. 
Hence (NW)(z) is 0(n~®~) for z on dU$ as n — > oo. For the inverse of 7VVF we 
get C(n e ) entries on dt/o- Hence the 0(n~i) term in the middle factor of the 
right hand side of (|4.87p turns into a 0(n _ 3) term: 

Po(z)N(z)- 1 =1+ 1 N(z)W(z)B ± 1 
712/(2)3 

x {'^ l ) ± {n^T{z))B ± W{z)- 1 N{z)- 1 + fn - ^) . (4.88) 

In the same way we find that the second term in the right-hand side of (|4.88|) is 
O(l) as n — > 00. We evaluate this term in more detail. 
By (|4~34|) we have 

N(z)W(z)B ± 1 = D^N{z)D{z)W{z)B ± l 

where by (|P5|) - ([P4"|) . P~6T]) . and (gTS]) , the last three matrices on the right-hand 
side are diagonal and their product satisfies for some constant c 7^ 0, 

D(z)W(z)B ± 1 = c Z P /3 {I + 0(z)) asz^O. 
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Furthermore the leading behavior of N(z) and its inverse as z — > follow from 
()4.49p and (|4.56|1 . Thus (|4.88[) reduces to (where we emphasize again the depen- 
dence on a), 

P Q (z)N(z)- 1 = 



niz3 $(z; a) 3 



x D^ 1 No{a){^i)±(n^T(z;a))N^(a)D 00 +0 (n~*) . (4.89) 
By (|437|) and (jOOl we have 

1 fV2i\ (7_ w i forlmzX), 

W0 ( a) = ^( /+O („ +1 ,,^gx|^ i ^ for[m2<o 

as a — > — 1. We plug this and the explicit formulas for ( v I / i)± (see Lemma l23|) into 
(|4.89[) . Then after some calculations we indeed obtain (|4.65p with Z n (z;a) given 
by gH. 

Since nr 2 (z;a) remains bounded as n — > oo if \z\ — n~ x / 2 and a — a n — 
— 1 + 0(n~ 1//2 ), we obtain that Z n (z; a) — 0(1). This proves the proposition. □ 

We note from (|4.84[) that Z n (z;a) is analytic in a punctured neighborhood 
of z = with a simple pole at z = 0. It also follows from (|4.84|) that 

Z n (z%; a)Z n {zi\ a) = for all zi, Z2 near 0. (4.90) 

This property will be important in the final transformations 



5. Final transformations 

We will do the final transformation S i->- R in two steps. First we will define Rq as 
the approximation error between S and the parametrices N, P_i,Po and Pi. The 
jump matrices of Rq will tend to the identity matrix as n — > oo on all parts of the 
jump contour except on 8Uq. 

Via a global transformation we finally define R, in such a way that it also 
has a jump that tends to the identity matrix on Uq. As a result of the steepest 
descent analysis, we then derive a global and uniform estimate for this function R. 

5.1. Transformation S n- Rq 

Define Ro(z) as 

S(z)N- 1 (z), z e C\ (SsUH^UCToUF!"), 

^ S(z)P -\z), Z€U \H S , 
_S(z)P^ 1 (z) 1 zGf/i\S s . 

Then i?o is defined and analytic on C minus the interval [a, 1], the lenses 
£^2, and the circles dU—i, OUq, and dU\. By comparing the jumps of S and 
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the parametrices we find that R has analytic continuation into each of the disks, 
and across the parts of the real intervals (a, 0) and (0, 1) outside of the disks. The 
singularities at a, and 1 are removable. For a and 1 this follows from the behavior 
of the Bessel parametrix given in [25]. For it requires a special check involving 
the behavior of * and ^P" 1 that we will not give here. The function Rq will then 
have jumps along the reduced contour E# shown in Figure [9] 

We choose clockwise orientation for the circles. The lips of the lenses are 
oriented from left to right, as before. Then Ro satisfies the following RH problem. 

• Rq is defined and analytic on C \ 

• i?o satisfies the jump relation i?o,+ = Rq.-Jr on ^R with 



N(z) J s (z)N- 1 (z), zeE fi \ {dU-i U dU U dUi), 

P^{z)N-\z), zeSU-!, f 

JR ° (Z)= <P (Z)N-Kz), zedUo, ( °' 2) 

P 1 (z)N-\z), zedUu 

• Ro(z) = I + (z^ 1 ) as z -> oo. 

Due to the matching conditions for the local parametrices (|4.6f [) and (|4.64[) 
we have 

Jro (z)=I + o(-) uniformly for z € dU-% U dUi (5.3) 



On the lips of the lenses the off-diagonal entries of the jumps of S involve 
the functions <fx,2, and they do not necessarily have negative real parts along the 
lenses S^ 2 - However, due to the estimates in Lemma 14.11 we may conclude that 
for a = a n = — 1 + ©(n" 1 / 2 ), we have for some constant C > 0, 

Re(pj(z;a n ) < ~Cn- 1/3 , for z G \ (dU-i U dU U dU\) 
Then by (j4TT71) we have J s = I + 0{e- CnV3 ), and so by JOJ) we also obtain 

Jr (z)=I + (e- Cn2/3 ^j as n oo (5.4) 
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with a possibly different constant C > 0. The O-term in (|5.4p holds uniformly for 
zeSji\ (dU-i U dU U dUi). 

Due to (|4.65p and ()5.2[) the jump matrix Jr on <9t/o is given by 

J Ro (z) = P (z)N(z)- 1 = I + Z n (z;a) + 0(n-^ for z e dU , (5.5) 

with Z n (z; a) given by (|4.84[) . and the O-term is valid under the assumption that 
a = a n = —1 + 0(n -1 / 2 ) as n — > oo, see Proposition 14.91 The jump matrix on dUo 
does not converge to the identity matrix as n —¥ oo. Therefore we need one more 
transformation. 

5.2. Transformation Rq h- R 

Let Zn\a) denote the residue of Z n (z; a) at the simple pole z = 0. Thus by (|4.84[) 



~(o), s r y , x 1 T(0;a)(nr 2 (0;a) + 9/3) 
>(a) = hm zZ n (z;a) = - — ^ 




(5.6) 



x D" 1 1 (V2i 1 -1) 5„ 



We define the new matrix-valued function _R as 

f 7 _^W 5 forzeC\(S fl UC7 ), 
^)^oWx « (0) (5.7) 

I f +Z w (g;o)- , forzeJ/o- 

Then R is defined and analytic in C \ £r, and satisfies 

• i? is analytic on C \ 

• R satisfies the jump relation R + = R-Jr on S/j with 

, M /(/ + Z B (,;a)-^)- 1 J JloW (j + ^) ,6^ 0) 

• = I + (z^ 1 ) as z ^ oo. 

Note that by (I4.84j) and (|5.6I) both Z n (z;a) and Z„ (a) are scalar multiples 

of the constant matrix D^ 1 (\/2i 1 — l) T (y/2i 1 — l) Doo whose square is 
zero. Therefore we have such relations as 

(l+Z n (z;a)-^^\ =I-Z n (z-a) + ^°± (5.9) 

and Z n (z; a) 2 = Z n {z; a)Z { °\a) = Z { ° ] (a)Z n {z; a) = (zi 0) (a)) 2 = 0. 

We also recall that for a — a n as in (|1.9I) we have by (|4.76l) that n 2 T(0; a n ) 
remains bounded as n — > oo, while f'(0,a n ) tends to y/2. Therefore by (15.61) 

Z^\a n ) =0{n- 1/2 ) asn^oo. (5.10) 
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This implies by (|5.2j) that the same estimates as we had in (|5.3|) . (|5.4[) for Jr 
remain valid for Jr on Er \ 9?7o- That is, 

Jl + O^- 1 ) forzedEZ-iUdC/!, 
Jfl(z) = \ J + O (e- c « 2/3 ) for zeE fi \ U dU U 9^). (5 " U) 



From (|4~65l) . ([O]) . (|5T8)) . and ([53]) we have when a = a n as in (fl~9 



Jr(z) = I- Z„(z;a„) + 



x (7 + Z„(,z;a„) +C (n b)) 



= / + Ofn-Bj, ze9C/o, (5.12) 
where we also used (|5.10p and the fact that Z n (z;a n ) remains bounded for \z\ — 

5.3. Conclusion of the steepest descent analysis 

We have now reached the goal of the steepest descent analysis. In the RH problem 
for R we have by (|5.1ip and (|5.12[) that all jump matrices Jr tend to the identity 
matrix as n — > oo. 

Then by standard methods, see e.g. [14] and also [10] for a situation with 
varying contours, we have that R{z) also tends to the identity matrix as n — > oo, 
at the following rate 

(n~ 6 \ 
I uniformly for z£C \ Sr. (5.13) 

This concludes the steepest descent analysis of the Riemann-Hilbert problem for 
Y. 



6. Proof of Theorem 11.41 

Recall that P n ,n{z; a) is the (1, 1) entry of Y, see (12.61) . which we write as 

P n , n (z;a) = (10 O)Y(z) . (6.1) 

The asymptotic formula ()1.10[) for P n>n will be derived from this by following the 
series of transformations Y n- T S m> R for z £ Uq. 

In the calculations that follow we assume that z is in the upper part of the 
lens around [0, 1]. The proof for other z in other regions is similar. We obtain from 
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([HP , 631, and jBH) 

p„ „( z;a ) = e 2n(9l(2)+92(z)) (1 0)S(z)[ j. (6.2) 

For z € J7 we have by (|4.78l) and (pTTI) 

S{z) = R (z)P (z) = R Q {z)E n (z ]a )y (nl/(z; a );n5r(z;a)) e nh ^W{z)-\ 

Inserting this into (I6.2[) and using the expressions (|4.66l) and (|4.67[) for A and W 
we obtain 

PnA z : a ) = z-^e 2n{9l{z)+32{z) e^ n{vi{z)+V2{z ^ (1 0)i? a) 

x * (n2/(z;a);n2T(z;a)) ^oj (6.3) 

where we used u>2(z) = z@(l — z)~>h 2 (z), see f|4. 12|) . The scalar prefactor in (|6.3p 
simplifies because of (|4.5I) and the result is 



PnA^a) = {-l) n e-^^ +l ^z-P (I 0)R o (z)E n (z;a) 

x * (fi5/(z;o);n5T(z;o)) ^oj (6.4) 

Note that by (|Q3|) , (|4T50]l , dHU, and 

( Z {a) (a)\ 
R (z)E n (z; a) = R(z) I + Z n (z; a) - - ^ ' D^N{z- a) 

x D(z)W(z)A (n%f(z; a);Jr(z; a)\ (6.5) 

Assuming that a — — 1 + C(n -1 / 2 ) and z — 0(n~ 3 / 2 ) as n — > oo, we have by (|4.56|) 
and (j4~57] 

i /^A 

^ ) = ^ 4 1 1 ^V + Ofa- 1 ' 2 ))- (6-6) 

Using this in the product (|6.5j) we see that the terms Z n (z; a) — Z " ^ are canceled 
by the leading order term in (|6.6I) because of the special form (j4.84j) of Z n (z; a). 

Since Z n (z;a) — Zn — is uniformly bounded, we find that the first four 
factors on the right-hand side of (|6.5|) combine to (we also use (|5.I3p ). 

(/ + 0(n- 1 / 6 ))^ I7 g^^-M \\ (-w 1 (6.7) 



Double scaling limit for modified Jacobi-Angelesco polynomials 



43 



The behavior of D(z) is determined by (|4~46|) . Then by (|4~67l) 



D(z)W(z) = Ce 



C1+C2 



23' 



32 



(a+ 7 )+/3 



- Z P /3 B+(I + 0{n- 1/2 )) 



where B+ is the diagonal matrix given by ()2.8|) . This factor also appears in the 
definition (|4.79p of A(z; r). Then the product of the last three factors on the right- 
hand side of (|6.5I) gives us 



Ce Cl+C2 e~e nT 



(z . af [Y 23^ / nif(z;a) 

V 2lT g±(a+ 7 )+/3 y Z 
// 3 A -1/3 



x S17 1 



1/3 



(I + Oin- 1 / 2 )) 



(nt/(z;o)) y 

Since (— u) 1 w 2 ) Q^ 1 = (0 l), we obtain by multiplying (|6.7|) and 



-Ro(z).E„(,Z;a) 



(^ e ci+C2g-gnT(^;a) 



2S" 



n2f(z;a) 



L-i. 

3 3 



6ir ■ 2 1 / 6 3 i(Q+ 7 )+/3 I z / 

x^ 1 I 1 1 (0 l)(I + 0(n-V 6 )) (6.9) 



Now let us now replace z and a by the scaled variables 



V2n 3 / 2 



On = -1 



rr/2- 



It then follows from ([67TO]) . and (|4J4)) - (|4J5)) that 

n5/(z„;a n ) = z + O f™~M , n^T(z n ;a n ) = t + (n~^ , 



(6.10) 



3.11) 



as tl — y oo , where the O-terms hold uniformly for z in a bounded set. Then we 
obtain from (16.91) 



i?o(Zn)-En(.2„;a„) = ^= i « 2 2 



6tt 



x D 



35(0+7)+/? 

1 j (0 1) (7 + 0(?i- 1 / 6 )) (6.12) 
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We insert (f6~T0]) . (jOTj) . and (|6~T2]) into (|6~4]). We also note that depends 
in an analytic way on a, so that by (|4.38p 

(D- 1 )^ =£> (oo;a n )- 1 = CT^l + ©(r^ 1 / 2 )). 

Then we obtain 

1 2 

PnA^a n ) = (-l)« e -^+^ ±± j-Jf 



37T 35(0+7)+^ 

xiz^O 1) *(z;r) ^oj (/ + ©(n^ 1 / 6 )) (6.13) 

The constants l± and I2 depend on a and therefore by (|6.10|) on t. One can 
show that 



e-Tl ,1+,J ) = (l?)" 6 "^™' 6 "^ ' + ,,, ' l4) 
Thus we proved (|1.10[) where C n is given by f|1.14[) and Q is 

Q(z) := iz~ p (0 l)*(z;r)^oj. (6.15) 

Looking at the definition (|2.14l) of ^(z;t) with < argz < ?, we find from 

flSZElJ) 

Q(z) - far" (e 2 ^q'l(z) + (6.16) 

Then by the integral representation (j2 . 13[) of qi and qi we find by easy contour 
deformation from (|6.16[) . that for Rez > 0, 

Q(z)=iz~ f} J t'P- 1 exp - ^ + ztj dt, (6.17) 

with To as in Figure [3] Making the change of variable zt h-> t, we arrive at the 
integral in ([l.lljl , which shows that Q from (16.151) indeed agrees with the Q defined 
by (II . lip in the theorem. 

This completes the proof of Theorem 11.41 

Remark 6.1. The property (14.90)) implies that (7 + Z n (z; a)) -1 = I — Z n (z;a) and 
so by (l4~65|) 

(I - Z n (z; a))P (z)N{z)- 1 =1 + ©(n^ 6 ) for z E dU . 

This suggests an alternative approach in which we redefine Po(z) with the extra 
factor I — Z n (z; a) on the left. Using this redefined Po in the transformation (|5.ip 
from S M- Rq, we would arrive at a RH problem for Rq that has all jumps close to 
the identity matrix, but which has a simple pole at z = 0. This can be included 
in the RH problem by providing a residue condition at z — 0. The simple form of 
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the residue matrix Zn (a) allows one to remove the pole again in a transformation 
Ro i— > R as in |18j . The resulting R is then exactly the same as before. 
We thank Percy Deift for this remark. 
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